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CN ; Observations of cosmic-ray electrons and positrons have been made 

r-- ! with a new baUoon-borne detector, HEAT (the "High-Energy Antimat- 

^ '. ter Telescope"), first flown in 1994 May from Fort Sumner, NM. We de- 

CL|' scribe the instrumental approach and the data analysis procedures, and 

O ; we present results from this flight. The measurement has provided a new 

to '. determination of the individual energy spectra of electrons and positrons 

' from 5 GeV to about 50 GeV, and of the combined "all-electron" inten- 

^ sity (e+ + e") up to ~ 100 GeV. The single power-law spectral indices 

H ■ for electrons and positrons are a = 3.09 ± 0.08 and 3.3 ± 0.2, respec- 
tively. We find that a contribution from primary sources to the positron 
intensity in this energy region, if it exists, must be quite small. 
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1. Introduction 

Electrons and positrons relatively 
rare component of the cosmic radiation, but 
are distinct from all other cosmic-ray parti- 
cles because of their low mass and the ab- 
sence of hadronic interactions with the con- 
stituents of the interstellar medium. Previous 
measurements of electrons (e"'" + e~) up to 
roughly 1000 GeV ( [Prince 1979| , [Nishimura et 



al. 198q , [Golden et al. 1984| , [Tang 1984[) have 
sh own that thoir intonoity amounto to about 



1% of Lhe flux of pruLuns around 10 GeV, but 
decreases more rapidly with energy (oc E~^'^) 
than the proton spectrum (oc E^'^-'^). Sepa- 
rate measurements of positrons and electrons 
have only been possible at much lower ener- 
gies and have indicated a "positron fraction" 
(e"*"/ (e^ + e~)) of a few percent in the region 



1-10 GeV ([banselow et al. 196c 


), 


Buflington 


3t al. 1975[, [Golden et al. 1996 


, 1 


Jarbiellini et 



d. 1996[) . Some measurements have indicated 



an increase above about 10 GeV ( [Agrinier 



t al. 1969[, [Miiller fc Tang 1987[, [Golden et 



d. 1987 , [Golden et al. 1994 ). These observa- 
tions have led to a number of conclusions, but 
have also left some key questions unanswered: 

1) The predominance of negative electrons 
can only be explained if electrons are accel- 
erated by primary sources. The alternate 
production mechanism, as secondary parti- 
cles from interstellar nuclear interactions of 
hadronic cosmic rays (mostly through the 
TT^ ^ /x^ ^ e^ decay), would yield neg- 
and positive oloctrons in about equal 



atii vo ana positive electrons m 
propui Lions, and thus can only account fui 
a small portion of the total e^ flux. It is 
worth pointing out that electrons are the only 
cosmic-ray component for which an extra- 
galactic contribution can be excluded with 
certainty: energy losses through Compton in- 
teractions with the 2.7K background radia- 



tion preclude their propagation through in- 
tergalactic distances. 

2) The steepness of the observed energy 
spectrum of electrons is usually explained as a 
consequence of radiative energy losses during 
propagation through the interstellar medium: 
inverse Compton scattering with photons and 
synchrotron radiation in the interstellar mag- 
netic flelds. Assuming that electrons are ac- 
celerated with the same energy spectrum at 
the source as nuclei (oc [Miiller et 



al. 199l| , [Swordy et al. 1993| ) the observed 



spectral slope is roughly consistent with the 
galactic containment time of ~ 10^ years 
of nuclei at GeV energies. However, there 
are several problems with this interpretation: 
First, the energy spectrum of electrons at the 
source is not known a priori, nor is it known 
that electrons and nuclei originate at the same 
acceleration sites. Second, it has been pointed 
out ( rrang 1984[) that the energy dependence 
(oc E~^'^) of the containment time observed 
for nuclei may be difficult to reconcile with 
the observed shape of the electron spectrum. 
Third, the "leaky box" assumption inherent 
in these explanations should not be applied to 
electrons, as it requires an unreasonably high 
density of electron sources in the galactic disk 
( [Gowsik fc Lee 1979[) . 

3) The small flux of positrons appears to 
be essentially consistent with a secondary ori- 
gin. If this is true, the production spec- 
trum of positrons can be calculated from 
the known flux of primary nuclei ( [Protheroi 



1982[) . Positrons would be produced contin- 



uously throughout the galactic disk, and the 
leaky box model could more easily be taken 
as an approximation for their propagation 
through the galaxy. In any case, a direct com- 
parison between the production spectrum and 
the observed positron spectrum would per- 
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mit a quantitative study of the propagation 
mechanism. However, it must first be ascer- 
tained that indeed all positrons are of sec- 
ondary origin. The previously observed in- 
crease in the positron fraction has led to much 
speculation about the appearance of primary 
positrons at high energy (e.g. [Harding fc Ka- 



ma jty 1987| , [Aharonian fc Atoyan 1991] , |Do 
gie |l fc Sharov^ggq , I'lyika 19891 , [rurner fc| 



Wi jlczek 19901 , [Kamionkowski fc Turner 1991| ) . 
Recent results ( Barwick et al. 1995 , Barwick] 
et ^1. 1997a| ) have not confirmed this increase, 
but have not yet shown conclusively that the 
positron flux is entirely free of primary con- 
tributions. 

It is clear from these considerations that 
many of the open questions can only be an- 
swered through measurements of electrons 
and positrons separately, and over as large 
an energy range as possible. This was the 
motivation for the construction of the High- 
Energy Antimatter Telescope (HEAT), an in- 
strument that utilizes a large superconduct- 
ing magnet spectrometer to separate positive 
and negative particles, and that incorporates 
powerful techniques to identify electrons and 
positrons and to reject hadronic background. 
A first balloon flight of HEAT was conducted 
in 1994. In the following, we shall describe 
this measurement and the data analysis tech- 
nique, and we shall present and discuss the 
results. 

2. Flight 

The flrst balloon flight of the HEAT e^ in- 
strument took place on 1994 May 3-5, from 
Fort Sumner, New Mexico, and data were col- 
lected at float altitude for about 29 hours. 
The payload reached a maximum altitude of 
36.5 km and drooped to a minimum of 33 km 
at night, as illustrated in Figure |I[ Over 



the course of the flight, the payload drifted 
between vertical geomagnetic cutoff rigidities 
of 4 GV and 4.5 GV, latitudes of 33.35° N 
and 35.3° N and longitudes of 100.0° W and 
104.3° W. The instrument was recovered un- 
damaged near Wellington, Texas. 

3. Instrument Description and Perfor- 
mance 



The expected intensity of cosmic-ray po- 
sitrons is quite low, of the order of 10~^ of 
the total cosmic-ray intensity at compara- 
ble energies. A successful measurement then 
necessitates a detector with large sensitive 
area to yield statistically signiflcant data, and 
with very high discrimination power against 
the overwhelming proton background. This 
is accomplished in the HEAT instrument, 
shown in Figure 0, through the combination 
of a superconducting magnet spectrometer 
(using a drift-tube hodoscope (DTH) tracking 
chamber) with particle identiflers employing 
a time-of-flight system (TOF), a transition- 
radiation detector (TRD), and an electro- 
magnetic calorimeter (EC). A thorough de- 
scription of the instrument appears elsewhere 
( [Barwick et al. 1997b| ), and the following just 
gives a brief summary of the apparatus and 
its performance. 



3.1. Time-of-FIight System 

The TOF measures the velocity of the par- 
ticle, distinguishes downward from upward- 
going (albedo) particles, and measures the 
magnitude of the particle's charge. The rejec- 
tion of upward-going particles is necessary as 
these would otherwise appear as particles of 
the opposite charge in the magnet spectrom- 
eter. Through the measurement of the mag- 
nitude of the charge, singly- charged particles 
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are identified and discriminated from lielium 
and lieavier nuclei. 

Tlie TOF consists of four scintillator slabs 
on top of the instrument, each with an active 
area of 100 cm x 25 cm, and of the top three 
scintillation counters of the EC (described be- 
low). Throughout the instrument, photomul- 
tiplier tubes (PMT) are used that are resis- 
tant to the magnetic fringe fields. 

The PMT signals are both charge- (ADC) 
and time- (TDC) digitized, and the parti- 
cle's charge is determined from the ADC val- 
ues measured with the top scintillator slabs, 
each of which has a PMT on either end. Fig- 
ure H shows the charge distribution obtained 
for flight data, with the criteria for selecting 
singly-charged particles shown as dashed lines 
(and summarized in Table |1|). The charge res- 
olution obtained is 0.11 charge units, which 
is sufficient to provide good rejection of He 
events. The TDC values, after corrections 
for the path length through the instrument, 
give the particle's velocity, with a resolution 
of 0.13 c (where c is the speed of light). Fig- 
ure ^ shows the velocity distribution obtained 
for flight data, and indicates the selection cri- 
teria used for the analysis. 

3.2. Transition-Radiation Detector 

The transition-radiation detector is used 
to distinguish electrons and positrons from 
hadrons. It is comprised of 6 modules, each 
consisting of a polyethylene-fiber radiator and 
a multiwire proportional chamber (MWPC). 
The MWPCs contain a gas mixture of xenon:- 
methane (70:30) and produce signals on both 
the anode wires and on cathode strips. The 
total charge deposited in the chamber is read 
from the cathode strips and pulse-height an- 
alyzed ("PHA" analysis), while the current 
signal is read from the anode wires in 25 nsec- 



wide time slices, which are digitized into three 
different threshold levels ("time-slice" analy- 
sis). 

Charged particles of high Lorentz factor 
(7 > 10^) produce transition-radiation (TR) 
x-rays (5-30 keV) in the radiator, which are 
subsequently detected by the MWPC. This x- 
ray signal is superimposed on the signal due 
to ionization energy loss of the parent parti- 
cle. In the energy range of 5-50 GeV, elec- 
trons produce a saturated TR signal, while 
protons and pions produce none, thereby per- 
mitting a discrimination between these spe- 
cies. 

3.2.1. PHA Analysis 

A maximum-likelihood technique is used to 
analyze the PHA signals ( [Cherry et al. 1974]) . 
To this end, first the probability distribu- 
tions for the pulse heights in the MWPCs 
must be determined, using clean populations 
of electrons and protons. The electron sam- 
ple is obtained by selecting events with nega- 
tive rigidity in the magnet spectrometer (de- 
scribed below), and a shower profile in the 
EC consistent with an electromagnetic shower 
(also described below). The proton sample 
is obtained by selecting events with positive 
rigidity and which have shower profiles in- 
consistent with an electromagnetic shower. 
The validity of this procedure has been ver- 
ified with accelerator calibrations at Fermi- 
lab. Gain differences between the TRD cham- 
bers are corrected for, as are temporal varia- 
tions (30%) in the chamber gains caused by 
temperature and pressure changes inside the 
instrument gondola, spatial non-uniformities 
(20%), and the relativistic rise of the proton 
ionization signal. The upper two panels of 
Figure ^ show the PHA probability distribu- 
tions Pe for electrons and Pp for protons. 
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For the subsequent data analysis, each e- 
vent, characterized by the MWPC pulse hei- 
ghts Xi {i =1,...,6), is compared with the prob- 
ability distributions, and a likelihood ratio 
LpHA is formed: 

L A ~ W ^'^^'^^^ 
i=l Ppi^i) 

Figure^ shows as a dashed line labeled "PHA" 
the electron and proton efficiencies obtained 
by applying a selection on the PHA likelihood 
ratio. 

3.2.2. Time-Slice Analysis 

The time-slice technique makes use of the 
fact that an x-ray photon absorbed in the 
MWPC generates a highly localized "ioniza- 
tion cluster" in the chamber, which then drifts 
to the anode wires and produces a "spike" in 
the time structure of the anode signal. In 
the analysis, again, a likelihood ratio L'rpg 
is constructed, using the appropriate single- 
chamber probability distributions Pg and P^. 
These distributions are based on the number 
of time slices above each threshold level, and 
the positions, heights, and number of clus- 
ters in the time-slice distribution, after the 
gain corrections described above are applied. 
The lower two panels of Figure ^ show the 
time-slice probability distributions obtained 
for electrons and protons. Figure ^ shows as a 
dotted line labeled "Time-Slice" the electron 
and proton efficiencies obtained by applying 
a selection on the time-slice likelihood ratio. 

The likelihood ratios Lp^A and L'j-g., are 
combined to form a total ratio: L Total = 
LpHA ■ L'rpg. Figure |^ shows the distribution 
of L Total together with the selection criterion 
used to identify electrons {L Total > 10^, or 
dashed line in the figure). With this crite- 
rion, electrons are retained with an efficiency 



eg = 88% and protons with an efficiency 
ep = 0.59%, corresponding to a proton re- 
jection power of 170. This is illustrated in 
Figure || ClS db solid line labeled "Total." 

3.3. Magnet Spectrometer 

The magnet spectrometer is used to mea- 
sure the rigidity R = ^ and charge sign of 
a traversing charged particle, and consists of 
a two-coil warm-bore superconducting mag- 
net and a precision tracking detector utiliz- 
ing drift tubes. The magnet produces an ap- 
proximately uniform field of central value 1 T. 
The fiducial volume within the magnet bore 
measures 50 cm x 50 cm x 61 cm. The 
tracking system consists of 479 drift tubes 
of 2.5 cm diameter, arranged in 26 rows; 18 
rows (of which 17 were operational during the 
flight) contain tubes parallel to the magnet 
axis, deflning the bending view, and 8 con- 
tain tubes perpendicular to the magnet axis, 
deflning the non-bending view. The drift gas 
mixture used is C02:hexane (96:4). 

The DTH determines the particle trajec- 
tory by measuring the drift times of the ion- 
ization tracks in each tube hit, from which, 
using a "time-to-space" function, the impact 
parameters (distance of closest approach to 
wire) are found. The trajectory of the particle 
in three dimensions is reconstructed from the 
impact parameters using a modifled version 
of the CERN program libraries' MOMENTM 
algorithm ( [Wind 1974| ) . Figures and |^ show 
the residuals from the tracks reconstructed for 
flight data, clS db distribution and as a function 
of radius, respectively. A signal is not used 
in the track flt if it lies more than one tube 
diameter away from the reconstructed trajec- 
tory, or if it has a residual that is worse than 
6 times the average residual for the track. 
Also removed are signals corresponding to an 
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impact parameter smaller than 2 mm, a re- 
gion where the tracking capabilities of the 
drift tubes are poor due to the nature of the 
electron-ion pair statistics and the increased 
drift speed near the wire. The single-tube res- 
olution achieved is about 70 fim. 

The performance of the magnet spectrom- 
eter is characterized by the maximum de- 
tectable rigidity (MDR). The MDR is defined 
as that rigidity where the RMS error in the 
sagitta measurement is equal to the sagitta of 
the track, and is computed on an event-by- 
event basis. The MDR distribution for elec- 
tron events, given in Figure |^, shows that a 
mean MDR of 170 GV is achieved. The MDR 
yields the relative RMS error in the rigid- 
ity: aR/R = R/MDR. Thus, an MDR of 
170 GV provides for a rigidity determination 
with an at least 3cr accuracy up to ~ 60 GV. 
In the analysis of the flight data, the selec- 
tions R/MDR < 0.25 and MDR > 60 are 
imposed to ensure meaningful results. 

3.4. Electromagnetic Calorimeter 

The electromagnetic calorimeter measures 
the energy of electrons and discriminates a- 
gainst hadrons. The EC consists of ten lead 
plates, each 0.9 radiation length thick and 
50 cm X 50 cm in area, and each followed 
by a plastic scintillator. 

Electrons deposit most of their energy in 
the EC, with well-understood shower profiles 
that provide a measure of their energy. Pro- 
tons and other hadrons, on the other hand, 
rarely interact in the EC, as it represents only 
about 0.3 proton interaction lengths, and if 
a nuclear interaction occurs, the longitudinal 
shower profile is usually quite different from 
that of an electron. Therefore, the EC also 
provides powerful discrimination against pro- 
tons. 



In order to reconstruct the electron pri- 
mary energy from the shower profile, Monte 
Carlo simulations have been carried out using 
the CERN program libraries' GEANT pack- 
age ( |Brun et al. 1994]) . The results of these 
simulations have been verified through accel- 
erator calibrations ( [Ibrbet et al. 1993| ). A 
covariance analysis is applied to the Monte 
Carlo events, providing the cross-correlations 
between the signals obtained in each of the 
layer pairs as a function of energy. Figure |T0| 
shows, for 10 GeV simulated electrons, the 
distribution of energies reconstructed by this 
covariance analysis. The correlation matri- 
ces generated by the covariance analysis of 
Monte Carlo events are then applied to flight 
data to determine the particle energies. The 
fractional energy resolution achieved is about 
10% for the energy range of interest. 

During the balloon flight, an event trig- 
ger is employed (see below) that requires a 
minimum energy deposition in the EC, and 
thereby removes most non-interacting pro- 
tons. In order to discriminate against the re- 
maining proton events with the EC, a mea- 
sure of the agreement between the observed 
and expected shower profile, generated by the 
covariance analysis, is determined for each 
event. In addition, the shower start depth t is 
derived from a separate profile analysis of the 
shower, and events with t > 0.89 radiation 
length are rejected. Finally, agreement be- 
tween the energy E measured by the EC, and 
the momentum p measured by the magnet 
spectrometer is demanded. Figure |TT| shows 
the electron and proton efficiencies obtained 
for flight data by imposing all of these re- 
quirements. The selection criteria chosen for 
the analysis are summarized in Table |I], and 
yield for the combination of shower counter 
and magnet spectrometer a proton rejection 
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power of 460 at an electron efficiency of 97%. 
Combined with that of the TRD (see above), 
the total rejection power of the instrument 
against protons is nearly 10^. 

3.5. Instrument Trigger 

The instrument trigger requires a through- 
going particle which deposits a signal larger 
than that of a 0.5 GeV electron in the EC. 
This excludes non-interacting protons. In ad- 
dition, a minimum number of hits in the DTH 
is required. A "fast" trigger is formed by 
the coincidence of signals in the top and bot- 
tom TOF scintillators, and of a signal above 
threshold in the sum of the lower seven EC 
layers. A "slow" confirming trigger is required 
within 1.5 /isec thereafter, based on the oc- 
currence of a three-fold majority of hits in 
the ffist, fifth, and last two rows of tubes in 
the bending plane. In order to obtain a sam- 
ple of non-interacting protons, the EC sum 
threshold requirement is removed for a small 
percentage (2%) of the events. 

4. Electron Selection 
4.1. Event Filtering 

The initial data analysis stage involves ex- 
amining the DTH hits for clean single-particle 
events. A rough estimate of the trajectory 
is made by performing a quadratic fit to 
the impact parameters in both the bending 
and non-bending views. Signals from tubes 
well outside this fit are discarded. Events in 
which tracks cannot be identified in the DTH 
through this procedure are rejected. Most 
events eliminated at this stage result from in- 
teractions of particles that penetrate the in- 
strument from the sides. 

Tracks which pass the initial DTH filter 
are then analyzed using the modified MO- 



MENTM algorithm, which determines the 
particle's rigidity from the measured track 
points. A parameter, generated by com- 
paring the measured track points with the fit, 
is used to reject events with large tracking er- 
rors. In addition, events are rejected based 
on the average residuals associated with the 
track, as well as on the number of tubes re- 
tained in the final fit in the bending and non- 
bending views. Finally, the MDR require- 



ments described in Section 2]3| are imposed. 
Table |^ summarizes all of the track cleanli- 
ness requirements. 

4.2. Template Fits 

Electron selection criteria are applied to 
events that satisfy the track cleanliness re- 
quirements. As discussed, these include selec- 
tions on the magnitude of the particle charge, 
the TRD total likelihood ratio, the EC covari- 
ance analysis parameter, the shower start 
depth, and reconstructed energy. In order to 
maximize the statistics of the results, the elec- 
tron selection criteria have been chosen rela- 
tively loosely, requiring that the number of 
background protons be less than 25% of the 
number of accepted positrons. Table |1| sum- 
marizes the selections used. 

For events meeting the electron selection 
criteria, distributions of the ratio of the mea- 
sured energy E determined by the EC, and 
the momentum p determined by the mag- 
net spectrometer, are generated. Figure ^ 
shows these E/p distributions for electrons 
and positrons in five top-of-the-atmosphere 
(ToA) energy intervals. (The top-of-the-at- 
mosphere energy EtoA is obtained from the 
measured energy E by correcting for atmo- 
spheric bremsstrahlung losses according to 
the procedure described in Section pl2| .) The 
differences in E/p between positrons and re- 
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si dual background protons is evident, with 
the positron distribution peaking at +1, and 
that for protons peaking at about +0.5. The 
tails at high \E/p\ values in the electron 
or positron distributions are due to brems- 
strahlung losses in the material above the 
DTH, which result in a reduced rigidity in 
the spectrometer while the entire energy (of 
the e^ plus bremsstrahlung photons) is usu- 
ally recorded in the EC. 

A fit to the E/p distributions for elec- 
trons and positrons is performed to obtain 
the number of electrons and positrons, and 
to estimate and subtract the residual pro- 
ton background. The E/p distribution for 
negative electrons, which can be taken as 
background free, is used to obtain a tem- 
plate of the expected E/p distribution for 
e^. The E/p distribution for background 
protons is obtained by inverting the TRD 
selection to select interacting protons (see 
Table |l|). Both template distributions are 
smoothed, and used to fit the E/p distribu- 
tion for positrons, shown as the solid curves in 



Figure |12| A Bayesian treatment is then used 
to estimate the electron, positron and proton 
counts, with prior probability distributions 
assumed to be fiat for the electron, positron 
and proton background counts, and with con- 
ditional distributions taken from the template 
fits. The electron and positron counts re- 
sulting from this procedure are listed in Ta- 
ble 0. The uncertainties given in the table 
are 16 and 84% Bayesian limits. For the en- 
ergy range 50-100 GeV, the ability of the 
magnetic spectrometer to determine the di- 
rection of bending is reduced as the uncer- 
tainty on the sagitta of the track is compara- 
ble to the sagitta itself. Therefore, in this en- 
ergy range, electron events are selected with- 
out the rigidity-related selection criteria of 



Table |T|, and electrons are not distinguished 
from positrons. 

As a test to verify the degree of potential 
residual contamination of the positron counts 
by interacting hadrons, the electron selection 
criteria were varied from very loose to very 
restrictive. In all but the loosest sets of selec- 
tion criteria, the electron and positron counts 
were reduced or augmented in the same pro- 
portion, indicative of little, if any, residual 
contamination. 

5. Absolute Energy Spectra 

The absolute differential energy spectra of 
cosmic-ray primary electrons or positrons are 
obtained from the raw electron or positron 
counts AA^ of Table S by calculating: 



Jpri 



AEeilAAt 



-Jsec{E) (1) 



with 



E 



Je/ E dE 
iS dE 



AE 



E-^'dE/E-'', 



where an atmospheric secondary (sec) com- 
ponent is subtracted to obtain the primary 
(pri) component. In equation (1), eQA is the 
instrumental acceptance. At is the live time, 
Ei and Ej are the lower and upper bounds of 
each ToA energy interval, respectively, E is 
the weighted average ToA energy in this in- 
terval, AE is the weighted ToA energy inter- 
val, froA is a correction factor related to the 
transformation of the energy scale to the top 
of the atmosphere, which will be described in 
Section f).2[ and a (~ 3.1 for primary parti- 
cles and ~ 3.2 for atmospheric secondaries) is 
the power-law spectral index. 
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The live time At of the detector is mea- 
sured in flight with a scaler which counts clock 
cycles only while the instrument is available 
for triggers and not busy processing an event 
following a trigger. Periods of instability of 
the instrument's performance are discarded, 
and a total live time of (18.04 ± 0.04) hours 
at float altitude is obtained. This live time 
includes losses due to transmission errors and 
data-unpacking errors from glitches in the 
electronics of individual subsystems. 

5.1. Monte Carlo Simulations 

Some of the quantities in equation (1) are 
determined with a full Monte Carlo (MC) 
simulation of the HEAT instrument, based 
on the GEANT software package, and includ- 
ing the actual detector configuration and the 
experimentally-determined fluctuations in de- 
tector response. For instance, the MC simula- 
tion is used to calculate the instrumental ac- 
ceptance for electrons, i.e., the product of the 
absolute efficiency e and the geometrical aper- 
ture QA. The MC-calculated efficiencies are 
compared to experimental quantities where 
possible, and in some cases renormalized to 
ensure agreement. For example, the TRD 
electron selection efficiency is determined by 
obtaining a set of electron events based on 
EC, TOF and DTH information, and mea- 
suring the fraction of such events that sat- 
isfy the TRD likelihood ratio requirement. In 
addition, a visual inspection of several hun- 
dred raw experimental and simulated elec- 
tron events revealed that about 65% of ex- 
perimental events are accepted by the anal- 
ysis, whereas 81% of simulated events are 
accepted. The instrumental acceptance de- 
termined by the MC is therefore corrected 
by a "scanning" efficiency factor of Cscan = 
0.65/0.81 = 0.80 ±0.08, the uncertainties be- 



ing an estimate of the spread of values due to 
independent scannings. 

Table |] lists, for each energy interval, the 
acceptances calculated with the MC, includ- 
ing the Cscan correction. The geometrical 
aperture, uncorrected for efficiencies, is about 
495 cm^ sr. This indicates an average elec- 
tron acceptance efficiency of about 37%. In 
Table |^, the acceptance uncertainties are ob- 
tained by adding in quadrature the uncertain- 
ties in the individual efficiency factors. The 
increased acceptance of the highest energy bin 
derives from an increased efficiency due to a 
different set of selection criteria used in this 
energy range, as indicated in Table |^. 

A comparison of experimental and simu- 
lated E/p distributions revealed a systematic 
energy calibration offset in the flight data. 
Specifically, the experimentally-determined 

\E/p\ peak for electrons is shifted upwards by 
about 14%, whereas the Monte Carlo sim- 
ulation predicts a shift of only 4% due to 
the bremsstrahlung energy losses. The 10% 
discrepancy appears to be due to a system- 
atic bias in the conversion of the ADC values 
of the EC into numbers of minimum-ionizing 
particles. This energy-scale shift is corrected 
for in Figure [T^ . 



The measured energy spectrum is shifted 
without a change in spectral index due to the 
finite, but roughly constant, energy resolution 
of the EC. The Monte Carlo simulation shows 
that of these events which have a true initial 
energy in a given energy bin, a fraction of 2 to 
6% are reconstructed by the analysis into the 
next higher-energy bin, while 22 to 32% are 
reconstructed into the adjacent lower-energy 
bin. This "spillover" effect reduces the recon- 
structed energy spectrum by about 10% in 
overall intensity as compared to the true pri- 
mary spectrum. The final results have been 
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corrected for this effect. Only events with 
ToA energies above 5 GeV, well above the ge- 
omagnetic rigidity cutoff at ~ 4.5 GV, are 
retained in the analysis. 

5.2. Atmospheric Corrections 

The electron spectrum measured at bal- 
loon altitudes is comprised of primary cosmic- 
ray electrons and atmospheric secondary elec- 
trons. The secondary component arises as a 
result of interactions of hadrons or primary 
electrons within the atmosphere, as well as 
reentrant albedo electrons. Reentrant albedo 
particles are significant only below the geo- 
magnetic cutoff energy. To determine the sec- 
ondary flux generated in the atmosphere, a 
Monte Carlo simulation utilizing the CERN 
libraries' FLUKA hadronic interaction algo- 
rithm (|Fass6 et al. 199^ ) is used. For this 
calculation, a primary proton spectrum with 
a power-law index of 2.74 is assumed 
t al. 199TD . Secondaries produced by pri 



eo 



mary electrons do not contribute significantly. 
To account for the contribution of heavy pri- 
mary nuclei, the simulated intensity of atmo- 
spheric secondaries is multiplied by 1.2 ( |Orth 
^ Bufhngton 19"76| ). The total intensity of 
10 GeV secondary electrons and positrons at 
6 g/cm^ is found to be 0.006 (GeV s sr m^)~^. 
This corresponds to overall corrections to the 
positron intensity of 20-30%, and to the elec- 
tron intensity of 1-2%, at the energies of in- 
terest. The absolute atmospheric intensities 
at an average atmospheric depth of 5.7 g/cm^ 
for electrons and positrons are also presented 
in Table |^, for each energy interval. 

As a test of the reliability of the MC cal- 
culation, atmospheric growth curves for elec- 
trons and positrons are measured during the 
flight. The positron fraction function 
of depth in the atmosphere is shown in Fig- 



ure 13 , for two separate energy intervals. The 
flight data are divided into four altitude in- 
tervals of 4.0, 4.5, 6.8, and 7.2 g/cm^ aver- 
age depth, respectively. The error bars rep- 
resent statistical fluctuations, dominated by 
the number of positrons in each energy inter- 
val. The positron fraction at the top of the 
atmosphere is obtained by linear extrapola- 
tion of the experimental growth curves. From 
the data on Figure |13|, one can determine the 
secondary-to-primary ratio {e~^ + e~)sec/ (e'*' + 
function of atmospheric depth. In 
Table ^ this ratio at a depth of 6 g/cm^ 
is given and compared to the corresponding 
results of the MC calculations. The uncer- 
tainties indicated are statistical. The MC 
simulation is consistent with the measured 
secondary-to-primary ratios within errors. 

The energy of a primary electron detected 
at a residual atmospheric depth of t radia- 
tion lengths is reduced from the energy at the 
top of the atmosphere due to bremsstrahlung 
energy losses. If the primary spectrum has 
the form of a power law the detected 

energy spectrum will, for a given depth t, re- 
tain the same power law but will be reduced 
in intensity by a factor This fac- 

tor is derived by taking the statistical fluc- 
tuations in the radiation loss mechanism into 
account ( [Rossi 1952| , [Schmidt 1972| ). For 



the HEAT flight described here, the atmo- 
spheric depth varies considerably throughout 
the flight. Therefore, we correct the energy 
for each detected electron or positron by a 
factor / = a*/"'"^^, where t is the current at- 
mospheric depth, and a = 3.1 is taken for the 
spectral index. Typically, this amounts to a 
~ 5-10% shift in the energy scale. Because 
equation (1) is used to calculate a differen- 
tial intensity, the energy scale shift necessi- 
tates a further correction factor JtoA = (/) = 
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(^^) = (1.086±0.020), where / is averaged 
over the flight. 

6. Results and Discussion 

The absolute differential intensities of elec- 
trons and positrons are obtained from equa- 
tion (1) using Table 0, which lists the num- 
ber of particles counted in each energy inter- 
val AA^, Table ^ which tabulates the accep- 
tai ^ces e QA, live time At and average ToA 



energy cuiiecLiun facLui JtoA) as well as E 
and /S.E (computed with a = 3.1), and Ta- 
ble ^ which gives the atmospheric contribu- 
tion. All of the statistical uncertainties asso- 
ciated with these quantities are also contained 
in these tables. The electron, positron, and 
all-electron differential intensities are listed 
in Table ^ The uncertainties are statistical, 
computed by adding in quadrature all of the 
contributing uncertainties. The results are 
plq tted in Figures pM| and p!5 | , scal e d by for 



cla iit)', Lugethei with ulhei meabUiemenLs uf 



electron and positron energy spectra 1 


Buffing- 


|ton et al. 1975, 


Nishimura et al. 1980 


, Golden 


3t al. 1984, Tang 1984|, Golden et al. 1994|). 



Besides the uncertainty in determining the 
instrumental efficiencies discussed in Section 
3TTI , there are two additional sources of sys- 
tematic errors. First, the necessity of apply- 
ing an energy-shift correction to ensure agree- 
ment with the momentum indicates a system- 
atic uncertainty in the absolute energy cal- 
ibration of a few percent, which translates 
into a ~10% uncertainty in the overall inten- 
sity normalization due to the a = 3.1 power 
law of the spectrum. Second, uncertainties 
in the normalization of the proton spectrum 
translate into a ~ 25% systematic uncertainty 
in simulated atmospheric secondary intensi- 
ties; this has a 0.5% effect on the electron 
and all-electron intensities, but results in 10% 



shifts in the positron intensity normalization. 
Thus, the estimated systematic uncertainty 
from these sources is 10% for the electron 
and all-electron intensities, and 14% for the 
positron intensity. 



In Figure |14|, we compare our results with 
previous measurements for which absolute in- 
tensities have been given for positrons and 
electrons separately. We find that the en- 
ergy spectrum of positrons reported by [Buh- 



ington et al. 1975| and [Golden et al. 199^ is 



generally consistent with our data, although 
our results show much improved statistical er- 
ror limits and permit a better estimate of the 
slope of the energy spectrum. The intensity of 
positrons of interstellar secondary origin can 
be calculated, provided the spectrum of par- 
ent protons is known, and the parameters of 
galactic propagation are specified. Such a cal- 
culation has been done by [Protheroe 1982[ , 
and, more recently, by Strong, Moskalenko 
fc Schonf elder 1997 . We show in Figure n 
the results of Strong et al. 1997 clS db dot- 
ted line, and notice the good agreement with 
our data. In Figure |15|, we compare the all- 
electron energy spectrum (e"*" + e~) derived 
from our data with previous results, also in- 
cluding measurements which did not employ 
a magnet spectrometer for charge separation. 
We arbitrarily chose to include only results 
published since 1980. The dotted line rep- 
resents a parametrization of the all-electron 
spectrum used by Strong et al. 1997. We 
find that all results exhibit similar spectral 
slopes over the energy range of concern, but 
that our overall intensity is lower than the av- 
erage of the previous investigations by about 
35%. There are two possible contributors to 
this systematic discrepancy. First, the abso- 
lute energy calibration for the experiments is 
likely to be uncertain by about 10%, leading 
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to a possible systematic uncertainty in the in- 
tensity of the order of 20%. Second, the as- 
sessment of the absolute detection efficiency 
of the instrument is notoriously difficult for all 
such detectors. Its determination does involve 
some intuitive judgment, although the more 
recent investigations, such as ours, benefit in 
this respect from more sophisticated Monte 
Carlo simulations than were previously avail- 
able. 

The individual spectral indices for posi- 
trons and electrons are a = 3.3 ± 0.2 and 
3.09±0.08, respectively, over the energy range 
5.0 to 50 GeV, as indicated in Table |^ and 
represented by solid lines in Figure |T1[ Thus, 
the positron spectrum appears to be slightly 
steeper than that of electrons. If all positrons 
are of interstellar secondary origin, one should 
expect that radiative energy losses eventu- 
ally lead to a power-law index that is larger 
by unity than that of the production spec- 
trum, which follows the ambient spectrum 
of the parent nuclei, i.e., a ~ 3.7 for nu- 
clear spectra of the form E~'^''^ . The spec- 
trum of electrons, on the other hand, should 
steepen less strongly, to a final slope with 



a : ^ 3.1 if electrons are mostly produced a~ 
the same primary sources as nuclei, and with 
the same source spectrum characterized by 
a ~ 2.1. However, these final spectral slopes 
may not be reached in the energy region be- 
low ~ 100 GeV. Thus, it would seem proper 
to fit our data not to a single power law but 
rather to a spectral form that refiects a tran- 
sition from the source spectrum (modified by 
solar modulation) to a spectrum that is fully 
steepened due to radiative energy losses. Such 
an analysis had previously been performed on 
the all-electron spectrum by Prince 1979] and 
[lang 198^ . Here we will not entertain this 
analysis for our data, as we soon should be 



able to improve the statistical quality of the 
results through the inclusion of data from an 
additional balloon fiight. Qualitatively, we 
may just conclude that the slightly steeper 
spectrum of positrons, as compared to that 
of electrons, is to be expected if positrons are 
predominantly secondary particles. This con- 
clusion must be confirmed by additional data 
and by an extension of the measurement over 
a wider energy range. 

As mentioned above, a comparison with 



the calculations of [Protheroe 1982| and Strong 
et al. 1997 indicates also that the absolute in- 
tensity of positrons is close to what may be 
expected if indeed all positrons are generated 
subsequent to nuclear interactions in the in- 
terstellar medium. More detail on the relative 
origin of electrons and positrons, and limits 
on the possible primary contribution to the 
positron fiux, can be obtained if the relative 
intensity of positrons, i.e., the positron frac- 
tion e"*"/ (e+ + e~) (which does not depend on 
the absolute efficiency of the detector), is in- 
vestigated as a function of energy. We have 
presented the results of such an investigation 
in previous papers ( |Barwick et al. 1995| , [Bar- 
wick et al. T997al) . 



The results described in this paper have 
been obtained with the first balloon fiight of 
the HEAT spectrometer. We believe that the 
high level of hadronic background rejection 
achieved with this instrument has led to data 
of high quality, limited only by the statisti- 
cal uncertainties of a single balloon fiight. A 
second balloon fiight which has been carried 
out successfully from Lynn Lake, Manitoba 
in 1995 August, covers the energy range from 
~ 1 GeV upwards. The analysis of the com- 
bined data set is presently in progress, and 
results will soon be reported. 
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Table 1 
Data selection criteria 



Cleanliness requirements 

DTH filter algorithm passed 

DTH rigidity reconstruction algorithm passed 

NB > 9 (Number of tubes in DTH bending view used in fit) 

NNB > 4 (Number of tubes in DTH non-bending view used in fit) 

< 10 (DTH rigidity goodness of fit) 

DEV-X < 0.080 cm (Average residual in X (non-bending) projection) 
DEV-Y < 0.014 cm (Average residual in Y (bending) projection) 
DEV_Z < 0.020 cm (Average residual in Z (vertical) projection) 
JBdl > 2.2 kG m (Integrated B-field over the track length)^ 
MDR > 60 GV^ 
MDR/ \R\ > 4^ 
\E/p\ < 3.0" 

0.70 < (5 < 1.65 (Velocity range, downward-going) 

\TOFJC - DrjC\ < 20 cm (Agreement between DTH track and TOF timing) 

Propagated DTH track traverses both TOF and EC 

Propagated DTH track traverses > 4 TRD chambers 

Electron selection criteria 

0.68 <Z < 1.45 (Charge) 

EC starting depth t < 0.89 radiation lengths 

< 2.6 (EC goodness of fit) 

logio [(^f/^ - ^LP^-^) • {Ll^ - cJlts)] > 3 (TRD hkelihood ratio) 



Interacting proton selection criteria 

Same as electron selection criteria except: 

logio [(^e7^ + <^L^»^) ■ {Ll' + ^LTs)] < -1 (TRD likelihood ratio) 



^This selection is not used in the 50-100 GeV energy range. 
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Table 2 
Raw numbers of electrons 



EtoA (GeV) # Electrons AA^e- # Positrons AA^e+ # All-electrons AA^e 



5.0 - 6.0 


1231 ±36 


107 ± 11 


1338 ± 37 


6.0 - 8.9 


1781 ±43 


161 ± 14 


1942 ± 45 


8.9 - 14.8 


918 ±31 




993 ± 33 


14.8 - 26.5 


340 ± 20 


18.51IJ 


359 ± 20 


26.5 - 50.0 


751^° 


6.11II 




50.0 - 100.0 






19.211:^ 
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Table 3 

Energy intervals and effective acceptances 



EtoA 


(GeV) 


E (GeV) 


/^E (GeV) 


e ^^4 (cm^ sr) 


5.0 


- 6.0 


5.45 


0.991 


181 ± 18 


6.0 


- 8.9 


7.16 


2.78 


193 ± 19 


8.9 


- 14.8 


11.1 


5.50 


194 ± 20 


14.8 


- 26.5 


18.9 


10.7 


189 ± 19 


26.5 


- 50.0 


34.5 


21.1 


160 ± 16 


50.0 


- 100.0 


66.4 


44.0 


207 ±21 



Note.— the live time is At = (18.043 ± 0.036) hr, 
and the average ToA energy correction factor is fxoA — 
(1.086 ±0.020) 
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Table 4 

Atmospheric secondary electron and positron intensities at 5.7 g/cm^ 



EtoA 


(GeV) 


j.Jm^ssrGeV) ^ 


j+im^ssrGeV) ^ 


5.0 


- 6.0 


(1.8 ±0.6) X 10-2 


(2.4 ±0.6) X 10-2 


6.0 


- 8.9 


(7.6 ±2.3) X 10-3 


(9.7 ±2.2) X 10-3 


8.9- 


- 14.8 


(2.1 ±0.6) X 10-3 


(2.5 ± 0.7) X 10-3 


14.8 


- 26.5 


(4 ± 1) X 10-^ 


(5.3 ±1.1) X 10-^ 


26.5 


- 50.0 


(7 ±2) X 10-5 


(9 ± 2) X 10-5 


50.0- 


- 100.0 


(1.1 ±0.7) X 10-5 


(1.2 ±0.7) X 10-5 
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Table 5 

Secondary-to-primary electron ratios at 6 g/cm^ 



Energy (GeV) 


MC sec./pri. 


HEAT sec./pri. 


4.5 - 6.0 


0.04 ± 0.01 


0.04 ± 0.03 


6.0 - 8.9 


0.04 ± 0.01 


0.04 ± 0.03 


8.9 - 14.8 


0.04 ± 0.01 


0.01 ± 0.03 


14.8 - 25.6 


0.05 ± 0.01 


0.01 ± 0.03 


25.6 - 50.0 


0.06 ± 0.02 
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Table 6 

Differential intensities of electrons, in (m^ssrCeV)"-^ 











E (GeV) 


jpriiE) 







5.45 1.13 ±0.12 

7.16 0.548 ±0.057 

11.1 0.141 ±0.016 

18.9 0.0278 ± 0.0033 

34.5 (3.64t°:i) X 10-3 
66.4 



0.076 ±0.016 
0.0405 ± 0.0070 
(9.21IJ) X 10-3 
(1.00t°:i) X 10-3 
X 10-4 



1.20 ±0.13 
0.589 ± 0.062 
0.151 ±0.017 

0.0288 ±0.0035 
(3.841°:^^) X 10-3 
(3.30lJi^) X 10-4 
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Table 7 

Power-law fits to the differential spectra 





a 


Jo (GeV-^ m-2 s-i sr) 


Electrons 


3.086 ±0.081 


227 ± 45 


Positrons 


3.31 ± 0.23 


24± 12 
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Fig. 1. — Altitude profile of the 1994 balloon flight 
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Fig. 2. — HEAT instrument schematic cross-section 
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Fig. 3. — TOF (5 and charge distributions 
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Fig. 4. — PHA and Time-Slice probability distributions for protons and electrons 
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Fig. 5. — TRD proton and electron efficiencies 
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Fig. 6. — TRD likelihood ratio L Total distribution for protons and electrons 
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Fig. 7. — Distribution of DTH residuals for bending-view tubes 
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Fig. 8. — DTH residuals as a function of tube radius for bending-view tubes 
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Fig. 9. — MDR distribution for electron events 
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Fig. 10. — EC energy distribution for 10 GeV MC events 
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Fig. 11. — EC proton and electron efficiencies 
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Fig. 12. — E/p distributions for electrons and positrons 
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Fig. 13. — Positron fraction as a function of atmospheric depth 
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Fig. 14. — Absolute differential energy spectra for electrons and positrons; the dotted line is a 
prediction from Strong et al. 1997, and the solid lines are power-law fits to the HEAT results 
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Fig. 15. — Absolute differential energy spectrum for all electrons (e"' 
parametrization from Strong et al. 1997 
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